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INTRODUCTION 


Atmospheric turbulence iii surface and low-altitude winds results 
in a number of important structural design problems to the engineer. 
This is especially true for the case of tall, thin structures such as 
smokestacks, towers and launch vehicles. This ground-wind loads 
problem on launch vehicles has been investigated almost totally using 
aeroelastic models in wind tunnels and applying techniques for 
predicting the full-scale load values. This approach was necessitated 
because of conflicting prelaunch operation requirements and by 
considerations for the actual flight vehicle. 

Because of characteristic differences between wind-tunnel flow 
and atmospheric surface winds, it was uncertain that the wind-tunnel 
results were accurately indicative of the full-scale loads. The wind 
tunnel presents a uniform velocity profile and a very low-turbulence 
environment for the model. Atmospheric winds near the ground present 
a nonuniform velocity profile^ due to shear flo\$ and turbulence values 
much greeter than those found in wind-tunnel flow. To study and 
evaluate the effects of turbulence and sheer flow on the response of 
launch vehicles, a research program utilizing a full-scale missile 
was initiated. 

A surplus Jupiter vehicle was erected at Wallops Island, Virginia, 
in an attempt to correlate wind-tunnel predicted loads with the 
measured full-scale values. The vehicle was instrumented to obtain 
base bending moments and tip deflection data on the freestanding 
launch vehicle while it was subjected to a range of atmospheric 
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surface winds. The full-scale ground -wind loads program is described 
and initial data are presented in reference 1. 

> photograph of the Jupi+ci- vehicle on the pad at the Wallops 
Island site is shown in Figure I. The missile is approximately 
60 feet in length and has a base diameter of 8.75 feet. Two wind 
sensors, capable of sensing orthogonal horizontal wind components arid 
following gust frequencies up to 5 cps, were installed on a mast near 
the Jupiter to help define the characteristics of the wind impinging 
on the missile. These two wind sensors were placed at heights of 
13 and 53 feet above the ground and can be seen in this figure. The 
instruments were located upstream of the vehicle for the prevailing 
wind direction at Wallops Island and far enough away from the vehicle 
to. avoid perturbed flew — on the order of five vehicle diameters. 

Since the wind sensors were this distance from the vehicle, the 
question arose as to a change in wind characteristics between 
measurement and impingement on the missile. One can extrapolate 
remote measurement along- the mean wind direction using Taylor's 
hypothesis which considers the turbulent velocities as a fixed field 
transported by and at the mean wind velocity. 

^ Under Taylor’s hypothesis, a space correlation function in the 
direction of the mean wind can be determined from the time correlation 
function using the transformation d = Ut. The spatial separation, d 
is measured in the direction of the mean wind, U. The nece/aary 
condition for the validity of this hypothesis is that the turbulent 
velocities have to be much smaller than the mean wind speed; i.e., 
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u/U « 1. It was felt that this hypothesis required experimental 
examination for these atmospheric ground wind studies. 

Along with this investigation of Taylor’s hypothesis, a 
significant amount of information on the statistical nature of low- 
level atmospheric turbulence was obtained. Some of this information 
is presented for its own merit and some is presented for comparison 
with values obtained in other investigations. 

Experience has shown that the local air velocities are continuous 
and random in nature and definable only in a statistical sense. The 
usefulness of expressing the properties of turbulence in statistical 
terms was first suggested by G. I. Taylor in 1921 (ref. 2) with three 
principal quantities of interest. They are (1) the relative frequency 
with which certain velocities occur, which is given by the probability 
distribution, (2) the frequency distribution of the energy contained 
in the wind, which is given by the power spectral density of the 
velocity fluctuations, and (3) the spatial correlation of th6 velocity 
fluctuations. The Gaussian distribution — a probability distribution 
identical to the Normal Law of Errors — generally gives an adequate 
description of the velocity distribution in atmospheric turbulence. 

On a log-log plot the spectrum of atmospheric turbulence shows no 
periodic motion and decreases linearly with frequency at a - * slope. 
These properties are examined for the wind data recorded at Wallops 
Island. The primary area of investigation, however, concerns the 
spatial correlation of the gust velocities as they are being 
transported along by the mean wind. 
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THE REVIEW OF LITERATURE 


The earliest at impts to provide a mathematical model to define 
turbulent flow assumed that the turbulent fluctuations consisted of 
discrete independent masses of fluid in random motion. The theories 
of this nature were developed mainly by Prandtl and resulted in the 
well-known mixing-length theory. It is now obvious that it was not 
realistic to consider discrete fluid particles which retain their 
identity over a certain distance. The discontinuous action implied 
by this theory was quite artificial, and the modern treatment of 
turbulence considers instead the continuous nature of the motion. 

The break from the theories that represented a discontinuous 
motion in turbulent flow was provided by G. I. Taylor in 1921 (ref. 2). 
The fundamental idea lay in the recognition that the velocity should 
be varying continuously with time along the path of the particle. 

Taylor again provided the next important advancement to the 
theory of the continuous nature of turbulent motion when in 1935 he 
considered the spatial structure of turbulence (ref. 3). This work 
introduced the correlation of velocities at two points as one of the 
describing quantities of turbulence. The statistical expression of 
this idea is provided by the cross -correlation function between the 
velocities at two points a given distance apart. 
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In this expression, homogeneous turbulence is implied in that the 
statistical properties R(d) and are taken to be independent 
of position (u^ = - u^). 

A time correlation, R(t), usually referred to as an auto- 
correlation function, may be defined in terms of the fluctuations 
measured at a point at instants separated by t. If the turbulence 
pattern passing over a point is unchanging and is being transported 
by the mean wind, it follows that 

R(t) = R(d) 

when d = U*t. This relationship was provided by Taylor in 1938 in 
the next important advance and is known as Taylor* s hypothesis 
(ref. 4). The necessary condition for this equivalence of time and 
space correlation functions through the transformation t = d/U is 
that the turbulence level in the flow be sufficiently low. The 
physical realization of this hypothesis is that, if the mean wind 
velocity is much greater than the turbulent components, the 
fluctuations at a point in space may be assumed to be the result of 
the whole turbulent field passing through that point at the velocity 
of the mean wind. The record of these fluctuations at a point will 
be nearly identical when measured along the axis of the mean wind. 

In wind-tunnel studies, under the necessary condition of 

u 2 /U tihc vailed ty ot T&y lor ' s hypotbosls h s B b66n duons 

The measurements of Favre, Gaviglio and Dumas (ref. 5) in a homogeneous 

flow indicate excellent agreement using the transformation t - d/U 

for a comparison between time and spatial correlation functions. 

<t 
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For atmospheric boundary layer flow the validity of the 
hypothesis is not quite so clear. The earliest observational data 
are available from an investigation by Giblett, et al., in 1932 at 
Cardington, England (ref. 6). These measurements were made at a 
height of 50 feet and provided some examples of simultaneous auto- 
correlations and space correlations from instruments approximately 
in line with the mean wind. Agreement between time and space 
correlation functions from these data is poor. 

A comparison of measurements taken from tower and airplane 
recordings at heights of 90 and 120 meters, respectively, conducted 
by Lappe, Davidson and Nctess indicated the hypothesis to be valid 
for the horizontal fluctuating components (ref. 7). 

An examination of Taylor's hypothesis at a height of two meters 
over smooth unobstructed grassland was reported by Panofsky, Cramer 
and Rao (ref. 8). The data considered for this study was recorded 
when the mean wind direction was within 10° of the line of wind 
sensors. Their work concluded that the hypothesis is valid at this 
height for horizontal separation distances up to 90 meters and 
intensity of turbulence levels as large as 0.26* 


THE EXPERIMENT 


Site Description 

Wallops Island, the launch site for Wallops Station, a facility 
of the National Aeronautics and Space Administration, is located on 
the Atlantic coast approximately 10 miles south of the Virginia-Maryland 
state line near Chincoteague, Virginia. It is comparatively flat, 
sparsely wooded, and sandy. Tho stationary tower used for measurements 
in this investigation was located about 300 feet from the shoreline and 
almost 500 feet from the nearest buildings. The surface surrounding 
this tower was flat and consisted of sandy soil with short grass as 
the only vegetation. For the high wind samples analyzed the wind came 
from the sea. A map of Wallops Island showing the shoreline, location 
of the stationary tower, and prevailing wind direction for the analyzed 
data is given in Figure 2. 


Apparatus 

The Wallops Island program obtained information on the statistical 
nature of atmospheric surface winds in general and the correlation of 
gust velocities along the mean wind vector in particular using the 
two-tower arrangement 3hown in Figure 3. A stationary tower on which 
four wind sensors were mounted at heights of 13 26, 39 and 53 feet 

above the ground was used in conjunction with a movable tower. The 
portable tower contained two wind sensors at heights comparable to the 
heights of the lower two sensors on the stationary tower (13 and 
26 feet). It was desirable to cover a greater range of heights, but 
to do so severely limited the capability of moving the smaller tower. 
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Fast-response drag-sphere anemometers are shown mounted at each 
level. 

Hie stationary tower is a vertical mast, circular in cross- 
section, approximately 75 feet in length, with a base diameter of 
20 inches and a diameter of 12 inches at the highest wind sensor 
level. The four wind sensors alined vertically on this tower are 
mounted on 6-foot horizontal booms which extend from the pole toward 
the northwest, as illustrated in Ficvr* 4(a). Each sensor was 
located three feet above its supporting boom. Nineteen markers for 
locking the portable tower were placed at 10° increments on an arc 
with a 40-foot radius from a marker directly beneath the wind sensors 
on the stationary tower. Each marker was identified as to its 
location relative to the marker below the wind sensors. For example, 
the marker located directly south of the reference marker was labeled 
180° and the one west 270°. A diagram of this arrangement is 3hown 
in Figure 4(b). 

The portable tower, shown in Figure 5, consisted of a 6-foot- 
square flatbed, two-wheel trailer on which a triangular structure 
was installed to support two wind sensors. The heights of the wind 
sensors on this tower were 13 and 26 feet above the surface. The 
trailer contained adjustable jacks and levels by which it was assured 
the tower was vertical and the instruments level. 

The drag sphere is a Langley-developed fast response wind sensor 
which is able to follow gust frequencies up to 5 cps. The instrument 
consists of a 7 -inch-diameter perforated hollow sphere mounted on a 
two -component force balance that senses orthogonal) horizontal 
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components of wind force. This instrument is described more fully in 
references 9 and 10. In addition to fast response, the drag sphere 
anemometer has several other important characteristics: (1) the 

natural frequency of the balance (60 cps) is high compared to the gust 
frequencies of interest; (2) the coefficient of drag is constant and 
independent of Reynolds number for both steady and unsteady flow in 
the Rjj range incurred for moderate atmospheric winds; (3) the drag 
force vector remains alined with the wind direction — a result of the 
perforations which stabilize the separation point for the flow about 
the sphere. Some results from measurements taken using this instrument 
are presented in references 10 and 11. 

A Bendix-Friez Aerovane wind transmitter was located approximately 
400 feet from the stationary tower and at a height of 65 feet above 
the ground. This instrument provided a continuous record of wind 
speed and direction on a Bendix-Friez strip chart recorder. 

The output signals of the drag spheres were amplified, FM 
multiplexed, and recorded on analog magnetic tape. The multiplexing 
enables one to record up to five separate signals onto one tape track. 
To recover the data the multiplexed signals are reproduced through 
discriminator units which convert the frequency modulated signals back 
to discrete voltages. 
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Procedure 

A typical data sampling interval was initiated by determining the 
prevailing wind direction from the continuous strip chart record of 
the Aerovane anemometer, placing the portable tower a known distance 
downstream of the stationary tower and recording data from all wind 
sensors on analog magnetic tape for a time interval of 15 minutes. 

The placement of the portable tower was accomplished using the markers 
for direction and a measuring tape for distance. Instantaneous gust 
velocities at two stations along the mean wind vector were thus 
simultaneously recorded. 
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Instrumentation 
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Figure 2.- Wallops Island test site. 
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Figure 3.- Two tower arrangement for low-level turbulence measurements. 













Figure 4*~ Test site description. 



DATA ANALYSIS 


The drag sphere anemometer measures orthogonal horizontal 
components of the dynamic pressures of the atmospheric winds. To 
determine wind direction from this instrument the orientation of the 
orthogonal axes has to be known. It is recommended that these 
sensing axes be alined north-south and east-west for field operation, 
and this was the orientation used at Wallops Island. 

Consider an instantaneous wind vector, U^, striking the drag 
sphere. is made up of a time-averaged mean value, U, and an 

instantaneous value of the fluctuations about the mean, u. . 

i 



This instantaneous wind vector strikes the drag sphere at some 
angle, 9^, relative to one of its sensing axes. 

A time history record of the 
output of the sensor gives a trace 


of orthogonal components of dynamic 
pressure, q n3 and q ew . Before 
each sampling interval a zero level 
signal and known calibrate signal 
were recorded to determine the 
scale for the recorded wind data 


N 



that followed. The procedure for S 

laboratory calibration and field 
use of the drag sphere anemometer 
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is given in the Appendix. These data were recorded on analog tape and 
had to be converted to digital values for computer analysis. 
Analog-to-digital transformation: 

analog tape sample rate = 10 samples/sec 
(10 samples/sec) (15 min) (60 sec/min) = 9000 points 
Thus, for and , i = 1 to 9000 

Dynamic pressure converted to velocity: 

At each sample time, t^, there are values of and 

i 

q from which the magnitude and direction of the vector quantity, 
ev. 

q^, is obtained using the following relationships: 



*i = N h 

The conversion from dynamic pressure to velocity is given by the 
following equations: 

(Oil = -m <]%]’ 

?i = 1 0,1 gi 

V = l°il 8ln 9 i 

u eu i = |»i| cos e i 
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The above 3teps are done for i = 1 -► 9000. The 15-minute time-average 
values of U and 6 are given by 


U - 




where 


8 = 


u 

ns 


= 1 
N 


N 


i=l 



N 

u =1 V 
ew N L-t 
* i=l 



and 


N = 9000 


The values now available are the mean wind speed, U, the angle, 
0, and 9000 instantaneous values of the north-south and east-west 
components of velocity. The objective is to determine the instan- 
taneous values of longitudinal (inwind) and lateral (crosswind) 
components of velocity. 

An axis rotation program using the following equations was used 
to give the longitudinal and lateral velocity components. 

Longitudinal, U: 

IT. = U cos 0 + U__ sin 0 
i ew. ns. 


Lateral, V: 



sin G + U 

ns 


cos 0 
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This essentially creates two new time history records in units of 

velocity (feet per sec) with the U component alined with, and the V 

component 90° from, the mean wind vector,. From these, 

N 

U = i ) U 
K L i 

i=l 


This value should be equal to U as given before axi3 rotation. 


N 



i=l 

The value of V should be zero when the correct 8 is used. 
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u- and v. are the instantaneous values of the fluctuations 
1 i 

about the mean in the longitudinal and lateral directions, respectively. 
The instantaneous values of these created time histories are then 
used for the statistical information obtained and presented in this 


thesis . 


RESULTS 


Seventeen sample intervals of 15 minutes duration each were 
recorded during ihe summer and fall of 1966 using the stationary and 
movable tower arrangement at Wallops Island. Seven of these sample 
intervals were chosen for extensive data analysis and are labeled 
data points one through seven for reference. The mean wind velocity 
as given by the Aerovane anemometer for each point is given in 
Figure 6 (a), while the wind direction and location of the portable 
tower relative to the stationary tower is shown in Figure 6 (b). The 
mean wind speeds for the seven points were between 37 (point 5 ) and 
48. 5 (point 2 ) fps, so a large range of mean wind speeds was not 
covered. These points were all recorded the 3 ame day and covered a 
total time interval of approximately four hours, during which time the 
weather did not change significantly. Therefore, the data were not 
influenced by changes in meteorological conditions. 

The ten other recorded samples were not statistically analyzed 
because the wind speeds encountered while recording them were less 
than 20 mph. The sensitivity of the drag sphere instrument is 
considerably reduced at low mean wind velocities. One reason for this 
reduced sensitivity at low wind speeds is that the drag sphere is, as 
its name implies, a drag-measuring instrument and its output is thus 
linearly related to the square of the velocity. Because of this 
reduced sensitivity at low wind speeds, it was decided to analyze 
extensively only t >5 samples that were recorded when the mean wind 
speed was above 20 mph. 
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The mean wind speed values as measured by the six drag sphere 
sensors for the seven data points are given in Table I. The letter P 
on the height value signifies the instruments on the portable tower. 

The maximum mean wind velocity measured at the 53-foot height was 
49.56 fps (point 2) while the minimum was 37.84 fps (point 5 ). The 
maximum value recorded at the lowest height was 39-87 fps (point 2) 
and the minimum was 28.19 fps (point 5 ). The values in parentheses in 
the 53-foot level column were obtained from the strip chart recording 
of the Aerovane anemometer output. 

The intensity of turbulence for each level is also given in 
Table I for the seven data points. This parameter is obtained as the 
ratio of the root-mean-square value (or standard deviation, o) of the 
longitudinal component fluctuations to the mean wind speed. It is a 
numerical measure of the amount of turbulence or gustiness contained 
in the flow. For the Wallops Island investigation, the maximum and 
minimum values of this uxameter were 0.183 and 0.084, respectively. 

The mean wind spea-is and intensities are presented as functions of 
height above ground in Figures 7(a) through 7(g). These plots indicate 
an increase in mean velocity and a decrease in intensity of turbulence 
with an increase in height. 

There are a number of analytical representations of the atmospheric 

wind profile available that are valid for certain meteorological 

conditions (refs. 12 and 13). In the lower part of the boundary layer 

of the atmosphere, an exponential approximation for the mean wind 

speed profile is frequently used. This power law approximation is 
given by 
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= v° /z i )n 

where 

U z = mean vrind velocity at height z 

U - mean wind velocity at height z 

z i 1 

n = nondimens ional quantity which depends on surface roughness 
and meteorological parameters 

The apparent reasons for the popularity of this power law representation 
are that it applies over a wide variety of meteorological conditions and 
that it is much easier to apply than the others that are more 
theoretically acceptable. 

Singer and Nagle (ref. 14) studied the variation of the exponent 

n with certain meteorological parameters such as wind speed, cloud 

amount, gustiness, lapse rate, time of day, and wind direction 

differences. They found that the values of n sh'/W a general decrease 

with increased wind speed up to about 20 ft/sec and then remain 

relatively constant. The Wallops Island d^oa were all above this 

value of wind speed. Also, since the Wallops Island data were 

recorded in a relatively short total elapsea time, it is felt that 

the other meteorological parameters did not change, and therefore the 

value of n should be constant. This is verified by the data. A 

curve with the exponent n = 0.17 best fitted the measured values of 

U for points 1, 2, 4 and 6, while n = 0.16 for points 3, 5 and 7 
z 

provided the best fit. This is in agreement with other Investigations 
of the boundary-layer profile over flat, unobstructed grassland (ref. 22). 


The mean wind speeds as measured from the stationary tower are 
plotted along with the mean wind speed values obtained from Aerovane 
anemometers located on a 250-foot meteorological tower on the island 
in Figure 8 for points one through four. These plots are interesting 
for two reasons: (1) the combined data show a continuous curve, and 

(2) the mean wind speed is still increasing with height at the 250-foot 
level. A photograph of this meteorological tower, which has Aerovanes 
located at 50-foot intervals up to 250 feet, is shown in Figure 9. 

Another form of presenting the velocity distribution in the outer 
region of the boundary layer is the velocity defect, U - U z . This is 
the difference between the free-stream velocity, U, and the velocity 
at a distance z above the surface. This velocity defect is 
normalized by u*, the friction velocity, and presented as a function 
of position in the boundary layer. This form is given in reference 15 
for turbulent flow in circular tubes and in reference 16 for turbulent 
flow over both smooth and rough plates. The form presented in 
reference 16 is reproduced along with the Wallops Island data in 
Figure 10. The distance above the surface is normalized by the 
boundary layer thickness. For z/5 < 0.15 the data presented in 
reference 16 are grouped around the logarithmic relation given by 
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For - > 0.15 the following empirical formula describes the 
experimental data 



Hinze states that apparently a value of z/S = 0.15 roughly indicates 
the boundary between the wall region and the outer region. 

The Wallops Island data were put in this form using a boundary 
layer thicicness of 253 feet and a free-stream mean velocity as the 
value measured at this height. The Wallops Island data for points 
one through four are closely grouped and are given by the dashed 
curve. The slope of this curve is in agreement with that of the other 
experimental data presented, but the values of the Wallops Island data 
are consistently lower by approximately 40 percent. 

The friction velocity, u ;;_ , mentioned previously i3 given by 

u " =v 'V p 


where 

= wall stress 

V 


p = density 
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This parameter is used in tha logarithmic profile approximation 
(ref. 1?) which is given by 

Z + T 

U = In 2 

z k z 0 

where 

k = von Karman constant (0.4) 
z Q = roughness length 

The roughness length, z Q , is a parameter characterizing the surface 
roughness. It can be described as the average height of the surface 
roughness projections and can v ary from 0.001 cm over ice (ref. 18) to 
several meters over cities. Using mean wind speed values from two 
/levels on the stationary tower and solving for u* and z Q , average 
values of u i; = 3 fps and z Q = 5 cm were determined for the seven 
data points. These values for the Wallops Island site are in 
agreement with the range of values for z q given in reference 18. 

A z^ = 5 cm places it in the long grass type of surface which is 
characteristic of the area on and around Wallops Island. 

The preceding information on profiles of wind and intensity and 
surface conditions was presented to describe the wind regime and 
environment from which the statistical parameters that follow W6re 


obtained. 


I 
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Any observed sequence cf values that cannot be described by an 
explicit mathematical relationship which gives the value of the 
quantity at some instant is called a random process. Therefore, 
since it cannot be described analytically, it must be studied and 
described in terms of probability statements and statistical averages. 
To analyze random data properly from sample time history records of 
finite length, it is necessary to establish three characteristics of 
the data. They are normality, stationarity, and randomness. All 
three characteristics were examined for the turbulence data recorded 
at Wallops Island and will be discussed at the appropriate time during 
the presentation of the results. 

Normality . - The probability distribution defines the relative 
frequency with which certain values of velocity occur. The Gaussian 
distribution, i.e., a probability distribution identical to the 
Normal Law of Errors, generally gives an adequate description of the 
velocity distribution in atmospheric turbulence. A comparison of the 
Gaussian form of distribution given by 

I _\-l -b/2a 

p(0 = ya v2tt J e 

and the distribution of the recorded data is given in Figure 11. A 
comparison is made for both the longitudinal, u, and lateral, v, 
component of velocity for a typical sample interval. The means, U 
and V, have been subtracted from the data so £ mean value of zero on 
these curves corresponds to the appropriate mean value of the signal. 
Since the mean has been subtracted from the data, the rms value is 
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equal to the standard deviation and the mean-square value is ual to 
the variance. 



u 


A visual examination of these curves shows close agreement between the 
measured distribution and the computed Gaussian distribution. 

The moments about the mean describe the shape of the distribution 
and provide another check on the normality of the data. The fourth 
central moment when nondimension&li zsd by the square of the second 
(o^) is called the kurtosis, or flatness factor, and has a value of 3 
for a Gaussian distribution. 

K = sV(o 2 ) 

The nondimensionalized third central moment is called the skewness, 
and if the distribution is symmetric about the mean, its value is zero. 

. , 2 3/2 

skewness = c 

The values of these moments are givt .n Table II and generally meet 
the requirements for the data to follow a normal distribution. 
Therefore, an assumption of a Gaussian distribution, normality of the 
data, for the u and v turbulent velocities appears to be valid. 
These values are in agreement with those obtained by Towns ond and 
reported by Batchelor in reference 19. Townsend’s work was with the 
u-eonponent turbulence generated by a wire-mesh grid in a wind tunnel. 
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If a process is shown to be Gaus.'ian, then the mean, th« mean 
square, and the auto covariance function completely determine its 
statistical properties. The autocovariance function defined by the 
equation 


C (t) = 


lim -3L 


T-fc-oo 


2T 


£(tk(t + 1 ) dt 


''-T 


describes the general dependence of the values of the data at one 
time, t, on the values at another time, t + t. The autocovariance, 
C^(t), is always a real-valued, even function with a maximum value at 
t = 0, and may be either positive or negative. The maximum value at 
i=0 is the moan-square value of the signal 


,2 = c (t = 0) 


Urn k 

T-*co 


^(t) 2 dt 


<-T 


The autocovariance is normalized by the mean-square value to give the 
autocorrelation function with maximum value of unity at v - 0. 


R r (-t) = C (t)/C (0) 

s s s 

The autocorrelation functions for the four levels on the 
stationary tower are presented in Figure 12(a) through 12(g). These 
correlation functions were obtained using increments of lag time, r, 
of 0.1 sec and with a maximum lag time of 12 seconds. These plots 
exhibit a dependence on height, with the breakdown in correlation 
occurring more rapidly at the lower levels. This is in accordance 
with the "long memory" and "short memory" regions of a wind-tunnel 
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boundary layer as concluded from experiments by Clauser and reported 
in reference 21. Clauser measured the wake generated by a rod placed 
at two distances from the wall, but still within the boundary layer. 
The disturbance decayed faster at the station nearer the wall and 
consequently a correlation value would be less at this station for all 
distances downstream of the initial disturbance. 

This height dependence of the autocorrelation function is also 
evident in a characteristic length of the atmospheric winds called the 
"scale of turbulence. " This scale of turbulence may be considered to 
be a rough measure of the largest distance that two points in a 
turbulent field may be separatsd before the correlation between gust 
velocities becomes zero. The autocorrelation function can be used to 
determine this characteristic length parameter* It is given by the 
product of the mean wind speed and the area under the autocorrelation 
function curve to where the curve crosses zero. This is given by 



where T Q = lag time at which R u (t) crosses zero. 

Since the autocorrelation functions presented do not cross zero 
for the maxi mum lag time (12 sec), "semiscales " were computed to 
determine the scale of turbulence. "Semiscalee, " as explained In 
reference 17, are defined as the lag distance at which the auto- 
correlation function dropped to a value of 0.6. Using the approximation 
of integral scales as twice the semiscales, the Integral scales so 
defined were computed by 
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L o.c 



R (i) d~ 
u 


where T_ , = lag time at which R (t) drops to a value of 0.6. The 
scales of turbulence are given in Table III and plotted in Figure 13. 

This characteristic length appears as a parameter in a 
mathematical description of turbulence (ref. 20) and reference 17 
reports evidence that for altitude below 1000 feet, its value is 
approximately equal to the altitude. For a linear approximation the 
scales computed from the Wallops Island data give a L = (4 x height) 
variation for the range of heights investigated. 

Homogeneity . - Figure 14 gives a comparison of the autocorrelation 
function at the same level as given by the instruments on both the 
stationary tower and the portable tower. The data are for points 5 
and 6 for which the instruments were separated horizontally by 70 and 
80 feet, respectively. The correlation functions for the stationary 
tower compare favorably with the corresponding function on the 
portable tower at the same height. This fact that the autocorrelation 
functions are essentially identical when measured at two points in 
space at the same height indicates homogeneity in the turbulence. 
Homogeneous turbulence is implied in that the statistical properties 
R (t) and u^ are taken to be independent of position. A complete 
test for horizontal homogeneity would involve measurements with the 
portable tower located at stations about the stationary tower other 
than those stations directly upstream or downstream of the stationary 
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tower. However, for the data presented it is shown that an assumption 
of homogeneous turbulence along the mean wind vector is valid for 
bounds on distance of at least 80 feet. 

Stationaritv . - Another characteristic of the data to be 
established is to determine whether it is stationary or not. For a 
signal to be stationary, its probability distribution and statistical 
averages are not functions of the time measured. A complete proof of 
stationarity would involve verification that all statistical properties 
are invariant with time. A practical test would be to divide a single 
sample and examine the mean, standard deviation, and autocorrelation 
functions for the parts. This has been done for a representative 
sample interval with the time history being divided into halves. The 
autocorrelation functions along with the corresponding mean and 
standard deviation values r re given in Figure 15(a) through 15(f). An 
examination of these functions ana values indicates that the data are 
reasonably stationary. 

Randomness . - The spectrum of the velocity fluctuations defines the 
frequency-wise distribution of the energy in the wind. The power 
spectral density function and the autocorrelation function are Fourier 
transform pairs for a stationary, random fluctuating quantity. These 
relationships are given by 

00 

PSD (f) = - I R (t) cos 2 nft dt 

U T7 I U 


0 


and the inverse 


r 


co 


K (t) = 


-'0 


PSD Cf) cos 2irft df 
u 


The autocorrelation function and the power spectral density function 
provide similar information in the time domain and the frequency 
domain, respectively. A test for randomness in the Wallops Island 
data would require a check for periodic or almost periodic components 
in the data. The most effective procedure for detecting periodic 
components in the data is an examination of the power spectral density 
function. A periodic component would be evidenced on such a plot by a 
noticeable spike at a given frequency. 

The power spectral density function i3 generally .1 malized by 
the mean-square value and this is the form presented in Figure 16(a) 
through 16(g). An examination of these curves shows no periodic or 
almost periodic components in the frequency range investigated — 

0 to 5 cps. Thus, the data are random in nature. The normalized 
power spectra as plotted in these figures show a decrease in energy 
content with frequency at a slope of - 5. This value of the slope is 
recognized as being valid for atmospheric turbulence spectra repre- 
sentations. The figures also exhibit no dependence on height for the 
spectra at the four levels on the stationary tower presented. 

Isotronv . - Isotropic turbulence is the condition for which certain 
statistical properties such as the rms values of the orthogonal 
components, u, v, w, of the gust velocities are unaffected by rotation 
of the reference axes, i.e. , they are equal and show no preference for 


direction. The rms values of the orthogonal components measured at 
Wallops Island, longitudinal u and lateral v, are given in Table IV. 
In general, the longitudinal component is greater than the corresponding 
lateral component with the differences in values being on the order of 
20 percent of their mean values. 

The primary area of investigation during this program was, as 
stated previously, to examine the characteristics of the turbulent 
velocities of the wind as they are transported by +he mean wind. A 
comparison of a fluctuating component at one station with its corres- 
ponding like component at another station along the mean wind vector 
is made by use of the cross -correlation function. 

Following the same procedure as in defining the autocorrelation 
function, we define the cross -covariance by 


C r (x) = lim 

**!» s 2 T->oo 



^i(tk 2 (t + 'O dt 


This relationship describes the general dependence of the values of 
one set of data, on the other, ^ le cross-covariance, 

C (t), is always a real-valued function which may either 

K V '2 

positive or negative. Also, CL r (t) does not necessarily have a 

^ 1 ’ ’2 

maximum value at x = 0 as does the autocorrelation function, or is 
it an even function. The cross -covariance is normalized by the 
product of the rms values of the individual signals to give the 


cross -correlation function 
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R (t) - C r (t 

‘■l' 4 ’1' "2 

with a maximum value of unity at some time lag. When the cross- 
correlation function does assume a value of unity at, for example, x, , 
it indicates that the two signals, £ and £ , are perfectly 
correlated —identical — with one lagging the other on the time scale 
by x, . When R (r) = 0, the two signals are uncorrelated, and 

1 h’ 4 

are statistically independent if this is true for all values of t. 

The cross-correlation functions were computed using the Wallops 

Island wind data for similar components measured at two stations along 

the mean wind direction. In other words, R (r) and R (t). 

u p u 2 v 1 ,v 2 

where 1 and 2 refer to the stationary tower and portable tower, 
respectively, were computed for both the 13-foot and the 26-foot 
heights. The convection time, i.e., the time it requires for the mean 
wind to convect the fluctuating velocities, u and v, the distance 
between the towers, is given by d/U. The maximum correlation value 
should occur when the lag time is equal to this convection time 
(t = d/U). 

The cross -correlation functions for the u component are 
presented in Figures 17(a) through 17(g). Figures 17(a) through 17(g) 
are in order of increasing separation distance between the two towers — 
from d = 20 feet to d = 100 feet. Each exhibits a lag time at 
which the correlation value peaks and this lag time to maxinpim 

, is compared to the convection time delay, d/U, in 

ID £UC 



correlation, t, 
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Figure 18. The values used to obtain this convection time delay are 
given in Table V. An examination of these figures indicates that the 
downstream signal does show a dependence on the signal measured at the 
upstream station which reaches a maximum at the appropriate convection 
time delay and this maximum dependence — or correlation — decreases 
with separation distance between the two measuring stations. 

Figures l?(a) through 17(g) also show that the correlation values are 
greater at the higher level instrument. This is to be expected since 
the intern ty of turbulence profiles show a decrease with height. This 
is also in agreement with Clauser's "long memory" and "short memory" 
regions of the boundary layer. 

A comparison of the time and spatial correlation values is given 
in Figure 19. The time correlation 

P- u (t) oc u^tju^t + t) 

is given by a representative autocorrelation function for the 
appropriate height. The space correlation is given by 

R (d) cc u.(t)u„(t) 
u l> u 2 x=0 ^ 2 

where and are recorded at two stations separated by a hori- 

zontal longitudinal distance d. The spatial correlation values, as 
represented by the circle symbols in this figure, were obtained from 
the cross-correlation curves (Fig. 17) at zero lag time. Although the 
space correlation values are in all cases less than the corresponding 
time correlation value, the data show good agreement between the two 
functions. These figures appear quite similar to data presented by 

p 

Ha . - - - - ' - 


Panofsky, Cramer and Rao in reference 3. From such data Panofsky, 
Cramer and Rao concluded that Taylor’s hypothesis, d = U", is valid 
within the limits of their investigation. 

From the Wallops Island data it appears that Taylor's hypothesis 
on the equivalence of time and spatial correlation functions does 
have some validity. The limits on the Wallops Island investigation 
are a horizontal separation distance of 100 feet, a height range up 
to 2b feet above the surface, and the intensity of turbulence in all 
sample intervals of less than 0.175. 

A more stringent interpretation of Taylor’s hypothesis would be 
to assume that the time history of a fluctuating velocity component 
does not change as it is being convected at the mean wind speed. A 
cross-correlation function between like components at two stations 
along the mean wind vectc would be identical to the autocorrelation 
function but shifted on the time lag scale by the appropriate 
convection time delay. In other words, the maximum value of the 
cross -correlation function would be unity and it would occur on the 
time lag scale at t = d/U. One would not expect this rigid inter- 
pretation of Taylor's hypothesis to be the :a3e in atmospheric boundary 
layer flow. The maximum cross-correlation values from Figure 17 are 
plotted in Figure 20 against the appropriate separation distance. If 
the time histories measured at the two stations were identical, the 
maximum cross -correlation value would be unity and independent of 
separation distance. This is given by the horizontal line at a 
correlation value of one. The experimental data should approach unity 
as the separation distance goes to zero. These curves do exhibit this 


characteristic with the values recorded at the 2 * -foot height greater 
in all cases than those recorded at the 13-foot height. Hie 
correlation values at a separation distance of ‘30 feet are low and 
apparently in error. The other data do decrease from unity as the 
separation distance increases. Therefore, this interpretation of 
Taylor's hypothesis, as expected, is not valid. 
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TABLE III.- SCALE OF TURBULENCE 


Point 

Height, faet 

13 

26 

39 

53 

13P 

26P 

1 

38. 28 

18.22 

144.89 

208. 67 

56.45 

112.34 

2 

47.60 

79.23 

119.47 

146. 52 

29.91 

42.29 

3 

51.50 

99. 37 

115.85 

155.07 

38.93 

51.87 

4 

57.12 

151.74 

254. 82 

240.60 

29.80 

58.16 

5 

44. 34 

74.13 

132.71 

227.51 

43. ;i 

54.52 

6 

47.16 

62.95 

122.33 

170.06 

43.02 

65. 96 

7 

56.00 

72.45 

103.19 

184. 14 

64.97 

117.02 























































TABLE IV.- RMS VALUES OF THE LONGITUDINAL AND LATERAL FLUCTUATING VELOCITIES 

(VALUES IN TABLE IN UNITS OF FPS) 


Height 

Data points 


2 

3 

4 

5 

6 

Tf 

i 

13' u 

1 

5.12 

5.19 

1 

4.52 

5.27 

5.59 

V 


4.01 

4.03 


3.42 

3.50 

4.23 

26' u 

wmm 

4.93 

4.77 

5.03 

5.10 

■Mi 

6.29 

V 

HUB 

4.24 

4.15 

3.43 

2.86 

H 

3.79 

39' u 

4.69 

4.34 

4.11 

4.29 

WBSM 

4.43 

4. 69 

V 

3. 66 

3.47 

3.33 

3.57 

1 

3 . 87 

4. 64 

5V u 

4.92 

4.45 

4.31 

4. 59 

3.53 

4.03 

4.32 

V 

3.31 

3.51 

3.75 

3.19 

2.67 

2.96 

3.66 

13' P u 

6.23 

3.72 

HfSSK 

5.34 

4. 68 

BH 

5.41 

V 

3.85 

3.79 

1 

3.88 

3.14 

in 

4.04 

26' P u 

5- 28 

4.64 


4.34 


I 

5.44 

V 

4.07 

4.31 

1 

5.08 

IPPPfPf 

1 

4.26 


f 





































TABLE V.. COMPARISON OF CONVECTION TIME WITH TDffi LAG TO MAXIMUM CORRELATION 


Point 

■si 


2 

d/0, 

sec 

Q 



13 

36.82 

37.42 

1.34 

1.6 

0.499 

1 

13P 

38.03 





(d = 50* ) 

26 

40.94 

40.82 

1.25 

1.2 

. 548 


26P 

40.71 






13 

39.87 

39.45 

1.01 

1.0 

.727 

2 

13P 

39.03 





(d = 40’ ) 

26 

43.77 

42.28 

.94 

.9 

.755 


26P 

40.79 






13 

38.45 

37.99 

.79 

.8 

.763 

3 

13P 

37.53 





(d=30‘) 

26 

42.64 

41. to 

.73 

.7 

.785 


26P 

39.52 






13 

34.90 

34.16 

.59 

• 6 

.739 

4. 

13P 

33.42 





(d » 20' ) 

26 

38.42 

36.50 

.55 

.6 

.796 


26P 

34. 58 






13 

28.19 

28.44 

2.46 

2.4 

.545 

5 

13P 

28.69 





(d =70’ ) 

26 

30.47 

29.81 

2.34 

2.3 

.681 


26P 

29.15 






13 

30.93 

31.30 

2.56 

2.9 

.475 

6 

13P 

31. 66 

. 




(d = 80' ) 

26 

34.75 

34.46 

2.32 

2.4 

.672 


26P 

34.16 






13 

31.19 

32.60 

3.06 

4.3 

.350 

7 

13P 

33.92 





(d = 100' ) 

26 

33.62 

• 34.25 

2.92 

2.7 

.565 


2tP 

34.88 





















































Mean wind speed, U, fps 
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O 1 PM 2PM 5 rM 4PM 5 PM 


Eastern standard time 
September 20, 1966 

(a) Mean wind speed variation during time interval in which 

points 1-7 were recorded. 
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Point 

G, deg 

d, ft 

1 

65 

50 

2 

65 

40 

3 

70 

30 

4 

65 

20 

5 

80 

70 

6 

80 

80 

7 

80 

100 


(b) Location of portable tower relative to stationary tower 

Figure 6.- Teat conditions. 




(a) Point 1. 

Figure 7.- Profiles of mean wind velocity and intensity of turbulence 

up to 53 feet above the ground. 
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(e) Point 5. 

Figure 7.- Continued. 










Figure 9.- Wallops Island 250-foot meteorological tower 
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(b) Point 2. 
Figure 12. - C< 
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(c) Point 3. 
ire 12.- Continued. 
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(d) Point 4 
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(g) Point 7. 
e 12.- Concluded. 
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(a) Point 5. 

Figure 14- - Comparison of autocorrelation function measured at the stationary- 

tower and at the portable tower. 
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(b) Point 6. 
Figure 14.- Concluded 
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(b) 26 feet. 
Figure 15. - Continued. 
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Autocorrelation 



(c) 39 feet. 
Figure 15.- Continued 
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Autocorrelation 



(e) 13 feet, portable tower. 
Figure 15.- Continued. 


Autocorrelation 
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(f) 26 feet, portable tower. 
Figure 15.- Concluded. 
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Frequency, cps 
(b) Point 2. 
Figure 16.- Continued. 
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(d) Point 4. 
Figure 16,- Continued. 
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(e) Point 5. 

Figure 16. - Continued. 






Normalized power spectra, 1/cps 
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(g) Point 7. 
Figure 16,- Concluded 


Cross -correlation 
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(a) d = 20 feet. 

Figure 17.- Cross -correlation function. 



(b) d = 30 feet. 
Figure 17.- Continued. 
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(c) d = 40 feet. 
Figure 17.- Continued 
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(d) d = 50 f 
Figure 17.- Cor 







Cross -correlation 
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(e) d = 70 feet. 
Figure 17.- Continued 
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(g) d = 1CX 
Figure 17.- C< 
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(a) Height = 13 feet. 

Figure 19.- Comparison of time and space correlation 
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(b) Height = 26 feet. 
Figure 19.- Concluded. 
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CONCLUSIONS 


from the data presented in this thesis, a number of conclusions 
can be made concerning the wind environment at Wallops Island, Virginia. 
The data were recorded using fast response drag sphere anemometers and 
covered a range of heights up to 53 feet above the ground. These 
conclusions concern the horizontal components of turbulence only 
since no attempt was made to measure the vertical component. The data 
presented resulted from seven data samples recorded during an elapsed 
time of approximately four hours. The mean wind speeds measured at a 
height of 65 feet were between 37 and 4S. 5 fps. The meteorological 
conditions did not change significantly during this time interval. 

The natural boundary layer studied was found to have the 
following characteristics: 

1. The mean wind speed profiles are given by the power law 
approximation 

u z = u z (z/ Zl r 

1 

with the exponent n having a value of 0.16 or 0.17. 

2. The intensity of turbulence decreases with height, and for 
the conditions of this investigation turbulence values were in the 
range of 0.084 to 0.183. 

3. The roughness length, z Q , is equal to 5 cm. This value is 
in agreement with other values obtained over flat, grass-covered 
surfaces. 
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4. The probability distribution of both the u-component and 
v-component appear to be Gaussian. 

5. The autocorrelation functions exhibit a dependence on height 
with the breakdown in correlation occurring more rapidly at the lower 
levels. 

6 . Hie "scale of turbulence" increases with height according to 
the linear approximation L = (4 x height) for the range of heights 
covered. 


7. The turbulence is homogeneous within bounds on distance of at 
least 80 feet. 

8. The turbulence is reasonably stationary. 

9. Atmospheric turbulence is a continuous random process. The 
normalized power spectra show a decrease in energy content with 
frequency at a slope of - ^ . 

10. Horizontal isotropy was not established. The rms value of 
the u-component is greater than that of the v-component on the order 
of 20 percent, as is the case in wind-tunnel boundary layers. 

From the cross-correlation data the following conclusions are 

made: 


1. The cross -correlation between like components at two stations 
along the mean wind vector indicates a dependence on height; i.e. , the 
correlation values are greater for the 26-foot height than for the 
13-foot height. 

2. The maximum cross-correlation value does occur at a time lag 
equal to the convection time delay, d/U. This maximum correlation 
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value decreases with separation distance between the two measuring 
stations. 

3 . The comparison between time and space correlations indicates 
the validity of Taylor's hypothesis, d = Ut, for separation distances 
up to 100 feet. 

4. A more stringent interpretation of Taylor’s hypothesis — 
assuming equivalence of the time histories measured at the two stations 
along the mean wind — is not valid. 
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APPENDIX 






DRAG SPHERE CALIBRATION FOR FIELD OPERATION 

The following steps are followed for each strain gage on the two- 
conponent balance to determine the gage sensitivity for use of the 
drag sphere instrument. 

1. Static calibration: 

Statically load the balance with incremental weights and 

record output voltage of obtained static calibration curve. 

+ i volts 

bs 


The curve should be linear with a slope, A^. A^ has units of 
lb/volt. 

2. Dynamic calibration: 

Place the instrument in a wind tunnel and record output 
voltages for various values of tunnel dynamic pressure. Multiply 
output voltages by to obtain the measured force in pounds and 
plot this force against the tunnel dynamic pressure. 

% 

•1 
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This plot should also be linear with a slope, A 2 . (The 
linearity of this curve indicates a constant drag coefficient, C^, 
in the Reynolds number range of calibration. ) 

F = (C D ) (drag area) (q) 
drag area = constant 

The product A^A 2 gives another constant, Ay with units of 
psf/volt. 

3. Calibrate resistor: 

With a known resistance, Rp in the circuit and no 
external load on the gage, the output voltage, V^, should be 
recorded: 


V^(volts) x A^ (psf/volt) = q^(psf) 


R^ is equivalent to a dynamic pressure of q^. 
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4. Field operation: 

With the drag sphere in place for field use and no wind on 

the instrument, insert the calibrate resistor, R , in line and 

record the voltage output, V . Since R is equivalent to q , it 

Fl 1 

follows that the sensitivity of the recordings will be 


A 4 = q i/V psf//volt 


FL^ is then removed from the circuit and the instrument is ready for 
field operation. 


